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Unitary operations acting on a quantum system must be robust against system- 
atic errors in control parameters for reliable quantum computing. Composite pulse 
technique in nuclear magnetic resonance (NMR) realises such a robust operation 
by employing a sequence of possibly poor quality pulses. In this article, we demon- 
strate that two kinds of composite pulses, one compensates for a pulse length error 
in a one-qubit system and the other compensates for a J-coupling error in a two- 
qubit system, have vanishing dynamical phase and thereby can be seen as geometric 
quantum gates, which implement unitary gates by the holonomy associated with 
dynamics of cyclic vectors defined in the text. 
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1. Introduction 

Nuclear magnetic resonance (NMR) has developed many techniques to control phys- 
ical systems and maintain their coherence [1, 2]. A composite pulse is one of such 
techniques, in which a sequence of pulses is employed to cancel out a systematic 
error inherent in the pulses [3] . A systematic error is an unwanted imperfection in 
control parameters, such as poor calibration, and should not be confused with a 
random noise. The composite 7r-pulse by Levitt and Freeman [4], developed with 
intuitive but convincing account of its robustness, opened up a new field of research. 
Now we have hundreds of composite pulses [5, 6] and dozens of methods to design 
them, such as iterative expansion [7], gradient ascent pulse engineering (GRAPE) 
[8, 9] and concatenation [10]. 

Recently, quantum information processing (QIP) [12, 13, 14, 15] has an influ- 
ence over the composite pulse design. Very accurate control of a quantum system 
is required for a successful quantum error correction, as shown in [14] for example. 
Any quantum algorithm can be simulated by quantum circuits composed of one- 
qubit unitary operations and the controlled-NOT (CNOT) operations. As a result, 
robustness is required for arbitrary one-qubit operations and CNOT operation. In 
contrast, operations with limited angles and phases have been required in conven- 
tional NMR manipulations. Numerous composite pulses have been proposed to date 
in the context of QIP [10, 16, 17, 18, 19, 20, 21, 22, 23]. 
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Geometric quantum computation [24, 25] has been proposed to attain reliable 
quantum control. In addition to the dynamical phase, cyclic evolution of a quantum 
system allows for various geometric phases [26, 27, 28, 29, 30], which are controllable 
and thereby can be utilised for unitary operations. We call a gate implemented with 
a geometric phase a geometric quantum gate (GQG) hereafter. Mathematically, 
a geometric phase is regarded as a holonomy associated with a closed path in 
a suitable base manifold associated with a cyclic evolution [31, 32, 33]. Random 
fluctuations along the integration path are expected to cancel out, leading to a 
quantum gate robust against random noise. Although there is numerical support 
for the robustness of GQGs [34] , this issue is still under debate [35] . 

In this article, we unite these two apparently different constructions of robust 
unitary operations. More precisely, we reveal that composite pulses robust against 
certain kinds of systematic errors are nothing but GQGs. This has been observed 
previously in one-qubit operations [36]. Now we elaborate and generalise this ob- 
servation to two-qubit operations, which arc indispensable for a universal set of 
quantum gates in QIP. Our work reveals that many composite pulses are geometric 
in nature and their robustness is attributed to the robustness of GQGs against 
certain errors. 

This article is organised as follows. Geometric phase, in particular Aharonov- 
Anandan phase and its application to implementation of a quantum gate are in- 
troduced in Sec. 2. We employ the perturbation theory as a guiding principle to 
design composite pulses and derive the robustness condition in Sec. 3. In Sec. 4, 
we present the main statement of this article, that is, existing composite pulses to 
suppress the pulse length error and the J-coupling error are GQGs. We will employ 
a group theoretical argument to present our statement in a unified manner. The 
assertion in Sec. 4 is exemplified in Sec. 5 and 6 by analysing various composite 
pulses from our viewpoint. Section 7 is devoted to conclusion and discussions. 



2. Geometric Quantum Gates 

Geometric phase, anticipated in many branches of physics and chemistry [26] , was 
formulated first by Berry in an adiabatic evolution of a quantum system. In [27], 
Berry considered a cyclic evolution of a quantum system whose Hamiltonian has 
time-dependent parameters, and pointed out that after the cyclic and adiabatic 
evolution, the system may acquire not only the dynamical phase factor, but also 
a geometric phase factor, which is given by a circuit integral in the parameter 
manifold. This integral is geometric, in the sense that it is independent of how fast 
the circuit is traversed. The Berry phase has been generalised in many ways. One of 
such generalizations is Wilczek-Zee holonomy: In the presence of n-fold degeneracy, 
the geometric phase factor can be replaced to an element of a unitary group U(n) , 
which is also independent of how fast the circuit is traversed [28] . 

Aharonov and Anandan showed in [29] that the geometric phase appears even 
in a non-adiabatic evolution. Consider an n-level system, whose normalised state 
vector at time t € [0,T] is given by \ip(t)) <G C". Dynamics of the system is charac- 
terised by the Schrodingcr equation 

i±m))=H(\( t ))m)), (2.i) 



Article submitted to Royal Society 



Geometric Aspects of Composite Pulses 



3 



where the Hamiltonian H(X(t)) is Hcrmitc and time-dependent through parameters 
X(t) = (Ai(t), . . . , Ajv(t))- Here we set h = 1. When the evolution is cyclic with a 
period T, i.e., 

\rb(T)) = e^(0)>, 7SK, (2.2) 

then the phase 7 the system acquires after the cyclic evolution includes geometric 
contribution 7 g , which is defined in terms of the dynamical phase 7a as 

7g = 7-7d, 7d = -/ dt(m\H(X(t))\Ht)). (2.3) 

JO 

This phase 7 g is called the Aharonov-Anandan phase. It is possible to interpret the 
Aharonov-Anandan phase in terms of geometric structure of the Hilbert space C n . 
See Appendix A for details. Also, for another expression of the Aharonov-Anandan 
phase, see, e.g., [29, 30, 33]. 

Applications of geometric phases are found in QIP. For example, Zanardi and 
Rasetti proposed to use the Wilczek-Zee holonomy to implement unitary gates [24]. 
It is also possible to implement unitary gates by using the Aharonov-Anandan phase 
[10, 25, 36, 37, 38]. To see this, let {|V'a)}i<a<n be the eigenvectors of a Hamiltonian 
H(X(0)) and suppose their dynamical evolution is cyclic, that is, 

\MT)) = U(T) \i> a ) , U(T) = Te-^o^K')) (2.4) 

and 

|Va(T)) = e i7 ° |Va> , 7 a eK, (2.5) 

where the time-ordered product is denoted by T. Equating Eqs. (2.5) and (2.4), we 
observe that \ip a ) is an eigenvector of U(T) with the eigenvalue e 17 , that is, 

U(T) |VO = e i7 ° |VO • (2-6) 
When there is no degeneracy, the spectral decomposition of U (T) is written as 

U(T) = {H + ---+ ^ l^n) (V-nl ■ (2.7) 

The phase 7" is decomposed as 7° = 7g + 7^ in terms of the dynamical phase 
defined as 

7d = - [ T dt(1>a{t)\H(\(t))\1> a (t)), \Mt)) = U(t)\*l> a ). (2.8) 
Ja 

A unitary operator U{T) is called a geometric quantum gate (GQG) if 7^ vanishes 
for all a. 

3. Perturbative Construction of Composite Pulses 

In actual situations in NMR, the dynamics is controlled by a sequential application 
of rf-pulses with constant field strength. Accordingly, the time interval [0,T] is di- 
vided into k intervals, in each of which the Hamiltonian is constant. More precisely, 
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we define the z-th temporal interval [ti-i,U}, where ti satisfies = to < t\ < ■ ■ • < 
ifc = T, and define a piecewise constant Hamiltonian, which takes the form H{X l ) 
in the i-th interval Here X 1 = (X\, . . . , X % N ) is a constant parameter vector 

while N is the number of control parameters. Then, the i-th rf-pulse gives rise to a 
unitary operator 

e~ iW \ W i = H(X i ) • (U - U-i), (3.1) 
and U(T) can be written as 

U(T) = e- iW " ■ ■ ■ e~ ml . (3.2) 

Now we wish to implement a 'target' unitary operator U as U = U(T). The target 
U should be implemented in a way robust against the error under consideration as 
much as possible. Hereafter we seek a condition for such robust implementation. 
We consider errors which cause displacement 

W i + 8W% (3.3) 

where SW 1 is a self-adjoint operator corresponding to the error. When 5W 1 is 
sufficiently small in the sense of the operator norm, we can use the perturbation 
theory and find 

e -KW+swn w e -iw* ( ln _ i6W fj . §W i := f dxe ixW *5W l e- ixW \ (3.4) 

Jo 

to the first order in 5W l . Here the identity operator on C" is denoted by 1„. The 
operator 6W^ is the error operator 8W 1 in the interaction picture. Then, the unitary 
operator U' implemented with the error SW % is given by 

U' = e-W+™*) . . . c -i(wi + swi) w jj (ln _ iAW) ; (3 5) 

where 

k 

AW = Y^V-^SWfV*- 1 , V 1 = e~ iWt ■ ■ ■ e~ iWl for i = 1,2,... k- 1,(3.6) 

i=l 

with V° — 1. Many, albeit not all, composite pulses satisfy the following robustness 
condition 

AW = 0, (3.7) 

which we can evaluate once we specify SW l . This condition guarantees the effect of 
the error vanishes to the first order in 5W l . 

Now we wish to address the relation between the robustness condition (3.7) and 
a classification of composite pulses common in the NMR community. There are 
two types, Type A and Type B, of composite pulses [5, 16]. The error tolerance 
is independent of the initial state vector for Type A composite pulses, whereas it 
is not the case for Type B composite pulses. In view of this, the composite pulses 
satisfying (3.7) are clearly of Type A. 
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4. Composite Pulses as Geometric Quantum Gates 

To see the geometric nature of Type A composite pulses, we follow the argument 
introduced in [36], which has been generalised to multi-qubit system in [10]. Suppose 
that the systematic error is proportional to W l : 

5W l = eW\ (4.1) 

As shown later, two kinds of systematic errors are of this form. The robustness 
condition (3.7) reads 

k 

AW = eY J V l ~^W i V 1 - 1 =0, (4.2) 
i=l 

where use has been made of the identity 8W{ = 5W l derived from Eq. (3.4) and 
Eq. (4.1). Taking the expectation value of AW with respect to \ip a ), we obtain 

k 

7d = E^ Q W = 0, 7 d Q M ~ - (Mi 1)\W* \Mi - 1)> . (4-3) 

i=l 

where \ip a (i)) := V 1 \tp a )- Hence, any composite pulse which is designed by the 
perturbation theory and compensates the error (4.1) is GQG. In what follows, we 
will show that composite pulses associated with two kinds of relevant systematic 
errors are GQGs. 

(a) Error on One-Qubit System 

We turn to a one-qubit system, whose Hilbert space is C 2 . An SU(2) operations 
we can implement with a single rf-pulse in NMR is limited to the form 

W l =e ini - a/2, (4.4) 

where rii = (cos</>,, sin^, 0) and er = (o~ x ,o~ y ,o~ z ) due to the apparatus limitation. 
Nevertheless, we can implement any SU(2) operation by combining at most three 
such pulses using the Euler angle decomposition [12, 15]. The displacement (3.3) 
under the error (4.1) is seen as 

0i->(l+e)0i. (4.5) 

This is a well-known systematic error called the pulse length error in the NMR 
community [5]. Hence, from the previous argument, we observe that any composite 
pulse compensating for the pulse length error is a GQG. 

(b) Error in Two-Qubit System 

For a two-qubit system, the relevant Hilbert space and the set of unitary opera- 
tions are C 2 ® 2 and SU(4), respectively In view of quantum information processing, 
the controlled-NOT (CNOT) operation 

£/cnotH0)(0|® Is + llHll®^ (4.6) 
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is important. f Here, \a) £ C 2 with a = 0, 1 is the eigenvector of a z with the 
eigenvalue (—1)°. The relevance of CNOT operation originates from the fact that 
any QIP can be implemented as a quantum circuit composed of onc-qubit unitary 
operations and CNOT operations [11, 12, 15]. 

By using the Cartan decomposition [15], CNOT operation can be rewritten as 
C/cnot = K\HK%, with 

H = e ™*<r*®<'* e i °'y°v®°v e i °"<'«®°* K U K 2 € SU(2) <g> SU(2). (4.7) 
Since a x ® a x is generated from a z ® a z through the following identity 

the Ising-typc Hamiltonian 

H = Ja z ® a z /4 (4.9) 

is essential to implement CNOT operations which is commonly realised in a weak 
coupling limit. Hereafter we shall be concerned with the J-coupling error defined 
by 

J (1 + e) J. (4.10) 

Several composite pulses robust against the J-coupling error have been proposed 
assuming that one-qubit operations are free from errors. These existing composite 
pulses [17, 18, 19, 20] are designed by making use of the following three generators 
only: 

X := cr z <g) a z , Y:=a z ®a x , Z:=l 2 ®<T y , (4.11) 

among the fifteen generators of SU(4). Evidently these operators satisfy su(2) al- 
gebra: 

[X/2, Y/2] = iZ/2, [Y/2, Z/2] = iX/2, [Z/2, X/2] = \Y/2. (4.12) 

Thus, we can construct an SU(2) subgroup by exponentiating the generators (4.11). 
Let us denote this subgroup by G. 
Now, let us put 

SU = J(U - ti_i)/2, 

W i = Cli {cos + sin^F) /2 = c-^ zl2 (fl l X/2) e^ z/2 . (4.13) 

Then, we observe that 

e -iw* = e -i* i z/2 e -in t x/a e i# < z/2 e q (4 14) 

Thus, for this W l , we observe the identification between f2j,X, Y, and Z and 
9i,a x ,(Jy, and <r z in Eq. (4.4), respectively. Since the J-coupling error (4.10) is 

f Precisely speaking, det (7cnot = — 1 and it is not an element of SU(4). Nevertheless, we can 
multiply this matrix by an unphysical phase e^/ 4 to make it an element of SU(4). Two quantum 
gates that differ by an overall phase will be identified hereafter. 
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equivalent to the pulse length error (4.5) under this identification, we can con- 
struct a "composite pulse", which is robust against the J-coupling error, if we 
merely replace 9i,a x ,a y , and o~ z by fli 7 X, Y, and Z, respectively. In fact, as stated 
before, such composite pulses based on the identification have been proposed in 
[17, 18, 19, 20]. One of these composite pulses shall be examined later. Composite 
pulses designed under this identification are GQGs, since this identification keeps 
the mathematical structure of the theory unchanged. 

Two remarks are in order. First, the definition of Z tells us that we can freely 
tune the parameter fa by changing the rf-field along the y-axis of the second qubit. 
Second, we can define the Bloch sphere for an orbit generated by G and \tp) <E C 2 ® 2 
if \ip) is an eigenvector of some element U E G. In other words, if there exists U E G 
such that 

C/|V) = e i7 |^}, (4.15) 

then the G-orbit G\fa) of \fa) is identified as the Bloch sphere S" 2 . This observation 
ensures that we can visualise the time evolution of a cyclic state associated with 
U G G as a trajectory in the Bloch sphere, as long as we use the composite pulses 
proposed so far. 



5. Examples of Geometric Composite Pulse 

In this section, we give several examples demonstrating our claim that two types 
of composite pulses introduced in the previous section arc GQGs. To this end, we 
shall evaluate the dynamical phase of several composite pulses and verify that the 
dynamical phase indeed vanishes in all cases. 

(a) One-Qubit System 
We parametrise our target U as 

U = exp (— \9n ■ tr/2) , n = (cos fa sin fa 0). (5-1) 

Then, from Sec. 2, a cyclic state \ip a ) associated with U is given as an eigenvector 
of U, that is, 

\iP a ) = \(-l) a n) , a = 0,1, (5.2) 
where |(— l) a n) is the eigenstatc of n ■ ex such that 

n.<T|(-l) o n> = (-l)°l(-l) n>- (5-3) 
We shall often use the following useful formula: 

(n\ m ■ tr\n) = n ■ m. (5-4) 

Note that the vector n is the Bloch vector for the state |n) and we have 

U\(-l) a n) =u; a \(-l) a n), u a = exp[(-l) a+1 i0/2]. (5.5) 
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All composite pulses, for which we evaluate the dynamical phases, are composed 
from k = 21 — 1 pulses, which satisfy the 'time-symmetric' condition 

W i = W h+i-i_ (5.6) 

Many implications of this condition are found in [5]. Now we address that this 
condition leads to 

72(0 =7d a (£ + !-*)■ (5-7) 
See Appendix B for the proof. Hence, the dynamical phase is rewritten as 

7d Q = 2 [ 7d a (l) + • • • + 7d a (Z - 1)] + 7d a (0 (5.8) 
for a composite pulse, which is made of k = 21 — 1 pulses. 

(i) 90°-180°-90° pulse 

The first composite pulse was proposed by Levitt and Freeman in 1979 based 
on a trajectory on the Bloch sphere [4]. This is a k = 3 symmetric composite pulse 
defined by 

e 1 = e 2 /2 = e 3 = %/2, ^ = ^ = 0, 4> 2 = tt/2. (5.9) 

We immediately find 

Wi = W 3 = (71-/4)4 a, W 2 = (tt/2)0 • a, (5.10) 

which leads to 

U = e - iWl e- iW2 e- iWl = -iuy. (5.11) 

Hence, we observe that the target is fixed to 9 = ir and (j> = ir/2 and there are no 
free parameters we may adjust. It follows from Eq. (5.11) that \ifj a ) = |(— l) a y). 
Let us proceed to the calculation of the dynamical phase. First, we have 

7d a (l) - -(vr/4) ((-lTy\ x ■ a \(-l) a y) = {-k / A)x -y = (5.12) 
from the formula (5.4). Next, we observe 

IV«(1)> = e-^ 1 m = l(-l) a y) = \(-l) a z) (5.13) 

to obtain 

7d a (2) = -(tt/2) ({-lfz\ y ■ a |(-1)^) = (-l) a+1 (7r/2)y ■ z = 0. (5.14) 
Summing up these, we reach 

7d a = 2 7d a (l) + 7d a (2)=0. (5.15) 
We can confirm Eq. (5.7) by further calculation. We find 

|Va(2)> = e- iw2 e~ iwl \t/, a ) = e~ iw2 \(-l) a z) = |(-l) a+1 i) , (5.16) 
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Figure 1. (Colour online) Excursion of the cyclic state \y) under the 90°-180°-90° pulse 
on the Bloch sphere. The green arrow is the Bloch vector of the cyclic state. The solid 
angle of the blue area enclosed by the trajectory of the Bloch vector is equal to n = 6, 
which also shows that this composite pulse is a GQG (See Appendix A). 

from which it follows that 

7d a (3) = -(tt/4) ((-l) a+1 z\ x ■ <x |(-l) Q+1 i) = (-1)«(tt/4)s ■ z = = 7d a (l). 

(5.17) 

The time-evolution of the cyclic states ends up with 

|Va(3)) = e~ iw3 |^(2)) = K-lfy) = \i/, a ) , (5.18) 
as expected. These results are summarised as 

— iW 1 — iW 2 — iW 3 

\±y) * ► \±z) 1 > \ T z) i > \±y) , 7d a ( 4 ) = 0. (5.19) 

See Fig. 1 for the graphical representation of this excursion. 

The lesson we learn from this composite pulse is that the converse of our state- 
ment is not always true: Not all GQGs for a spin- 1/2 system are Type A composite 
pulses robust against the pulse length error. Indeed, this pulse is of Type B since 
AW 7^ 0. This has been overlooked in [36]. 

(ii) SCROFULOUS 

SCROFULOUS is a k = 3 time-symmetric composite pulse constructed by 
Cummins, Llewellyn and Jones [21]. This composite pulse was designed by using 
perturbation theory and quaternion algebra. Given a target (5.1), SCROFULOUS 
takes the form 

6»i = 6» 3 = arcsinc[2 cos(0/2) /tt], 9 2 = n 

4>i = 03 = (ft + arccos[— it cos 6\/{26\ sin(0/2))], 

<t>2 = 4>i - arccos[-7r/(26»i)], (5.20) 
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Figure 2. (Colour online) Excursion of the cyclic state \x) of the SCROFULOUS for a 
target — n, (f> = on the Bloch sphere. The green arrow is the Bloch vector of the cyclic 
state. The state e~ lW \x) pauses during the application of the pulse W 2 , since it is an 
eigenstate of the pulse W 2 . The solid angle of the blue area is equal to 8; this composite 
pulse is a GQG. 



where sine a; = sinx/x. Note that SCROFULOUS implements any one-qubit uni- 
tary operator of the form (5.1). 

Let us evaluate the dynamical phase. We set <j> = for simplicity, while extension 
to an arbitrary <j) is straightforward. First, we have 

7g(l) = -e x {{-lfn\n v c7/2\{-l) a n) = (-l) a+1 ^n-m/2 = (-l) a+1 (0i/2) oos&. 

(5.21) 

Next, observe that 

V 1 W 1 = ■— [cos 2 (6<i/2)n2 + sin 2 (6>i/2)m - sin0i(m x n 2 )] ■ cr, (5.22) 
where 

m = 2(n 2 ■ TOi)tzi — n 2 . (5.23) 

Since we have 

n 2 ■ n = cos (/) 2 = — 7rcos</>i/(2#i) + smfa sin[arccos(— 7r/(2#i))], 

m ■ n = cos(20i — fa) = — tt cos fa/(26\) — sin^i sin[arccos(— 7r/(2#i))], 

(m x n 2 ) • n = 0, (5.24) 

we observe 

7d a (2) = -(^ a |U 1 t^ 2 V- 1 |^) 

= (-l) a (7r/2){(7r/(26'i))cos0i + cos6>i sin[arccos(— k/(26x))} sin <£i}. 

(5.25) 
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Using sin(arccosa;) = \/l — x 2 and 

sin0! = y/l- (7r/(40i)) 2 /sin(0/2), (5.26) 
we immediately derive 

7d a -2 7d Q (l)+7d a (2) 

= (-lfie^l - (tt/(20 1 )) 2 ]cos0i + (tt/2)cos0i sin[arccos(— k/(26 1 ))] sin^} 
= (-l) a [l - (tt/(29 1 )) 2 }[9 1 cos0i + 7rcos0i/(2sin(0/2))] 
= 0. (5.27) 
Hence SCROFULOUS is a GQG. The trajectory of the cyclic state is given in Fig. 2. 

(iii) Broad Band 1 (BB1) 

Now we turn to the BB1, which was proposed by Wimperis [39]. For brevity's 
sake, we treat a k = 5 time-symmetric variant of the BB1 sequence. We call this 
variant time-symmetric BB1. The BB1 pulse sequence is useful for the implementa- 
tion of QIP, since it compensates for the pulse length error up to the second order in 
perturbative expansion [16]. There are two techniques to generalise the BB1 pulse 
sequence [22]. Using these techniques, we can design a composite pulse sequence, 
which compensates for the pulse length error up to an arbitrary higher order in 
perturbative expansion. 

For a target (5.1) with angles 9 and 4>, the time-symmetric BB1 consists of 

9 1 = 9 5 = 9/2, 9 2 = 3 /2 = 9 4 = tt, 

01 = 05 = 0, 02 = 04 — + K, 03 = 30 + K, (5.28) 

with 

k = arccos[-0/(47r)]. (5.29) 

Let us evaluate the dynamical phase associated with the time-symmetric BB1. First, 
we note from U = e~ 2lW that 

V 1 |VO = e~ iwl |VO = ±^T a |<M ■ (5.30) 

Then, we have 

7d a (l) = - W 1 |V ) = {-l) a+1 9 ini ■ n/2 = (-l) a+1 9/4. (5.31) 
Next we find from 9i = tt and <p2 = + K that 

7^(2) = - (Va| V 1 V 2 !/ 1 |VO = - (Va| W 2 |Va) = (-l) a+1 7m 2 • n/2 = {-!)% 
and 

e _iM/2 |VO = l(-l) a n') ; n' = (cos(0 + 2k), sin(0 + 2k), 0). (5.33) 
This leads to 

7S(3) = - (V' a | U 2 V 3 U 2 |Va) = (-l) a+1 7r n 3 • n! = (-l) a 0/4. (5.34) 
By adding individual dynamical phases, we finally obtain 

7 | = 2 7d Q (l) + 2 7d X2) +71(3) - (-l) a+1 (9/2 - 9/4 - 9/4) = 0. (5.35) 
This result confirms that the time-symmetric BB1 is also a GQG. 



(5.32) 
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(iv) Knill's sequence 

Knill's sequence [40, 41] is a k = 5 time-symmetric composite pulse. This se- 
quence implements the target U given by 

6 = 7r, n = (cos(a — 7r/6), sin(a — 7r/6), 0) (5.36) 

where a is a free parameter. The sequence is defined by 

9i = 7T (1 < i < 5), 01 = 5 = Q + 7r/6, 02 = 04 = Oi, <p 3 = a + 7r/2. 

(5.37) 

This sequence is used in experiments to maintain the coherence of nitrogen- vacancy 
centres in diamond [40] and to decouple a system from the environment [41]. Note 
that this sequence is robust against not only the pulse length error, but also the 
off-resonance error [41]. 

Let us calculate the dynamical phase. First, we have 

7^(1) = - (V> | W 1 |Va) = (-ir +1 nn v n/2 = {-l) a+1 / 2) cos(7r/3) = (-l) a+1 7r/4. 

(5.38) 

We find V 1 |V>a) = |(-l) a w') with 

n' = (cos(a + tt/2), sin(a + tt/2), 0). (5.39) 
Then, by the similar argument as that used for the first step, we have 

7d(2) = 0. (5.40) 
Further, we observe V 2 \i/) a ) = |(— l) a n") with 

n" = (cos(a - tt/2), sin(a - tt/2), 0). (5.41) 

Then, we have 

7d a (3) = (-1)V2- (5.42) 
We find, by adding individual dynamical phases, 

7d Q = 2 7 g(l) + 2 7d Q (2) + 7d Q (3) = (-1) Q+1 (V2 + - tt/2) = 0. (5.43) 

This example shows that the composite pulses robust against several systematic 
errors are also GQGs, if they compensates for at least the pulse length error. Thus, 
by construction, the composite pulses proposed in [10, 23], which are simultaneously 
robust against the above two errors, are also GQGs. 

(6) Two-Qubit System 
Since our interest lies in the CNOT operation, we choose the target 

U = e-' lQX/2 , (5.44) 

which is the entangling part in the CNOT gate. The cyclic state |^ a ) is an eigenstatc 
of X in Eq. (4.11). In the binary notation a = 2p + q where p,q £ {0, 1}, we find 

IVO = \P) ® \q) • (5.45) 
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Jones designed a composite pulse sequence for a two-qubit system from a one- 
qubit composite pulse sequence [17], by employing the isomorphism among the 
generators given in Sec. 4 (b). Let us introduce a notation (O)^ = exp[— ifl(cos(j)X + 
sin0y)/2] and set the target to (tt/2)q in this notation. Jones' sequence is given by 

(7r/4)o(7r)«(27r) 3K (7r) re (7r/4)o, k = arccos(-l/8). (5.46) 

Since the isomorphism maps X, Y, and Z to the Pauli matrices a x , <r y , and cr 2 , 
respectively, Jones' sequence is a two-qubit analogue of the BB1 sequence: the 
combination of the first and last pulses is the target pulse (9 = tt/2, <j> = 0) and the 
others are the same as the BB1 sequence (5.28). Similarly, the composite pulses in 
[18, 19, 20] are the two-qubit counterparts of those in [22]. 

Evaluation of the dynamical phase is easy if we make use of the isomorphism 
already mentioned. Since X is mapped to o~ x , the cyclic vector \p) <g> \q) should be 
sent to |(— l) p+q x) : which is also an eigenvector of the target U — cxp(— iira x /A). 
Thus the dynamical phase of Jones' sequence is transferred to that of the BB1 
sequence, which leads to 

7d = 0, (5.47) 

showing the sequence has vanishing dynamical phase. One can also achieve the 
same result by direct calculation without employing the isomorphism. 

6. Two Composite z-Rotations 

In NMR, rotations around the z-axis must be implemented by a sequence of pulses, 
since the rf-pulses (4.4) have the restriction rii _L z. Thus, it is of interest to 
investigate whether the sequences are geometric. 

First, we consider the following k — 3 sequence to realise a target U = or l6a ~l' 1 : 

9 1 = 6 3 = tt/2, 9 2 = 9, 0! = -03 = tt/2, 2 = 0. (6.1) 

The cyclic states are \ip a ) = \ (~^) a z) = \ a )- Let us calculate the dynamical phase. 
The first one is 

7§(1) = -(tt/4) m y ■ a Ha) = (-l) Q+1 (^/4)y ■ £ = 0. (6.2) 
We find V 1 \tp a ) = \(-l) a x), which leads to 

7d a (2) = -(9/2) ((-iyx\ x ■ a \{-lYx) = {-l) a+1 9/2. (6.3) 
Furthermore, we obtain V 2 \tp a ) — cxp[(— l) a+1 i0/2] |(— \) a x). Thus, we observe 

7 $(3) - (tt/4) ({-l) a x\ y ■ a \{-l) a x) = (-l) a (7r/4)y ■ x = 0. (6.4) 
We conclude 

7d ' = (-l) a+1 9/2 0. (6.5) 

Hence the pulse sequence (6.1) is not a GQG. Note that this sequence is not robust 
against the pulse length error, that is, AW ^ 0, which is exactly the contraposition 
of our claim. 
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(a) (b) 

Figure 3. (Colour online) Excursions of the cyclic states on the Bloch sphere, (a) The 
trajectory of the cyclic state \z) under the pulse sequence (6.1). Note that the trajectory 
fails to close, which shows that this sequence is dynamical, (b) The trajectory of the cyclic 
state |— z) under the pulse sequence (6.6) for 6 = ir. The solid angle subtended by the 
trajectory of the Bloch vector is it = 8, which shows the geometric nature of the sequence. 



Second, we investigate a k — 2 pulse for U = Q iea ~/ 2 \ 

9i = 9 2 = 7T, fa = 0, fa = 6/2. 
The cyclic states are the same as those of the previous sequence. We have 

7 g(l) = -(tt/2) (il> a \x-a H'a) = (-l) a+1 (n/2)x ■ z = 0. 
By the same way, we compute 

7d a (2) = 0, 

which clearly shows 

7d = 0. 



(6.6) 



(6.7) 



(6.9) 



Hence the pulse sequence (6.6) is a GQG. This pulse is not robust against the pulse 
length error. Indeed, we may check 



AW 



1 + cos — I a x — sin — a v 



7^0 



(6.10) 



by direct calculation. This also tells us that not all GQGs are robust against the 
pulse length error. The difference of these two composite z-rotations are visualised 
in Fig. 3. 



7. Conclusion and Discussions 

In this article, we uncovered the relation between GQGs and the composite pulses 
robust against certain kinds of systematic errors. For the error (4.1), proportional 
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to the Hamiltonian times the operation time, the compensation of the error auto- 
matically leads to vanishing dynamical phase. Thus, a non-trivial operation by a 
composite pulse robust against such an error is a GQG. 

We pointed out that there are two kinds of errors assuming the form (4.1). One 
is the pulse length error and the other is the J-coupling error. This implies that 
the composite pulses robust against these errors are GQGs. This observation was 
illustrated and confirmed by directly showing that the dynamical phase vanishes 
for several typical composite pulses: 90°-180°-90°, SCROFULOUS, BB1, KnilPs 
sequence for the pulse length error and Jones' pulse sequence for the J-coupling 
error. The two most common composite z-rotations were also examined. 

Our work has shown that we can construct a universal gate set composed of 
GQGs simply by using the composite pulses. This suggests that NMR is quite a 
useful test bench of geometric quantum computation. In view of this, further study 
of composite pulses, e.g. [37], is desirable for deeper understanding of the geometric 
quantum computation. 

We would like to thank Jonathan A. Jones and Yukihiro Ota for valuable discussions 
and Dieter Suter for drawing our attention to KnilPs sequence. This work is supported 
by 'Open Research Center' Project for Private Universities: matching fund subsidy from 
MEXT (Ministry of Education, Culture, Sports, Science and Technology), Japan. YK and 
MN would like to thank partial supports of Grants-in-Aid for Scientific Research from the 
JSPS (Grant No. 23540470). 

Appendix A. Geometry of Aharonov-Anandan Phase 

In this appendix, we outline the relevant aspects of the Aharonov-Anandan phase in 
the context of the present article. The geometric nature of the Aharonov-Anandan 
phase is derived from that of the fibre bundle structure associated with the Hilbert 
space. See [13, 32] for technical details. 

Consider the Hilbert space C n . In quantum mechanics, we are exclusively con- 
cerned with the set of normalised vectors in C™ . The set of normalised vectors form 
the (2n — l)-dimensional sphere S* 2 ™ -1 C C™. Moreover, we need to identify vectors 
that differ by an overall phase; two normalised states \if>) and e 17 represent the 
identical physical state for any 7 e I, The manifold obtained from S 271 ^ 1 under 
this identification is called the complex projective space; 

cp"- 1 ~ s* 2 "-7u(i), 

where U(l) is the set of overall phases. 

For n = 2, we obtain CP 1 = S 2 , which is nothing but the Bloch sphere. Ac- 
cordingly, S 3 is identified with a U(l)-bundle over S 2 (the Hopf fibration). More 
generally, S 271-1 is a U(l)-bundle over the base manifold CP™ -1 . A point in CP™ -1 
represents a physical state and its phase freedom is represented by the fibre U(l). 
The identification naturally introduces the projection tt : S 271 ^ 1 — > CP" -1 . Fixing 
the phase is equivalent to taking a point in the fibre (See Fig. 4). 

It should be noted that CP™ -1 has a natural metric called the Fubini-Study 
metric. Given a metric in the base manifold, we can construct a connection in the 
base manifold. This defines the horizontal lift of a given curve in the base manifold 
to the fibre bundle S' 2 ™" 1 . Now suppose that there is a closed loop in the base 
manifold. If one carries a point on a fibre over p G CP™ -1 along the horizontal lift 
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U(l) 




Figure 4. Schematic diagram of Aharonov-Anandan phase. The set of the normalised states 
forms S 2 "- 1 C C n and this subset S 2 "" 1 can be seen as a U(l)-bundle over the projective 
Hilbert space CP™ -1 . Given a closed path in the base manifold CP™ -1 , the horizontal 
lift of the path is naturally defined by a connection in the U(l)-bundle. The holonomy 
associated with the horizontal lift is given as e 17g €U(1), which can be seen as a global 
phase difference accumulated through the parallel transport along the horizontal lift on 
the path in CP™ -1 . 



of the loop, the point comes back to a point in the same fibre, which is not necessary 
the initial point. This U(l) phase factor obtained after traversing a loop is called the 
holonomy associated with the loop and the horizontal lift. The Aharonov-Anandan 
phase is nothing but this U(l) phase factor, which is geometric in the sense that 
it depends only on the loop in the base manifold and the connection of the Un- 
bundle but not on how fast the loop is traversed. 

We note that the twice the Aharonov-Anandan phase is the solid angle at the 
origin subtended by the trajectory of a state vector on the Bloch sphere (CP 1 , in 
this case) during a 1-qubit operation. 

Appendix B. Proof of Eq. (5.7) 

In this appendix we prove Eq. (5.7). For this purpose, we first note the identity 

a z W i a z = -W i (Bl) 

for W % of Eq. (4.4), because rii _L z. Multiplying — i and exponentiating Eq. (B 1), 
we find 

a z e~ iw 'a z =e iw \ (B2) 

Then, we obtain 

V k ~ i = e iw * ■ ■ ■ e iWl JJ = a z V i a z U (B3) 
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for a time-symmetric composite pulse. It then follows that 

\1>a(k-i)) = V^^a) 

= * t V*<r t U\il>a) fromEq. (B 3) 
= uj a a z V l (J z \ij) a ) from Eq. (5.5) 

= LOaO-zV 1 IV'ael) 

= u a a z \ip a<sl {i)) , (B4) 

where we denote the sum modulo two by ©. Therefore, using the condition Tr W z = 
and the completeness relation with respect to {\i/> a (i))} a =i.2, we observe that 

7§(fc + l-») = -(Mk-i)\W k+1 - l \Mk-*)) 

= - (V'oeiWkzWVz |V>offii(»)) from Eq. (B4) and |w a | 2 = 1 

= (VoffllCtJlW^IVaffllCt)) fromEq. (B 1) 

= ^[^(la-^aWX^-WI)] 

= - (^a(i)l |<Aa(i)) from TrVU 1 = 

= -<tf a (i- l)|W*|Va(i-l)) 

= 7d a W, (B5) 
which proves Eq. (5.7). 
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